Disordered jammed packings of frictionless spheres 
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At low volume fraction, disordered arrangements of frictionless spheres are found in un-jammed 
states unable to support applied stresses, while at high volume fraction they are found in jammed 
states with mechanical strength. Here we show, focusing on the hard sphere zero pressure limit, that 
the transition between un-jammed and jammed states does not occur at a single value of the volume 
fraction, but in a whole volume fraction range. This result is obtained via the direct numerical 
construction of disordered jammed states with a volume fraction varying between two limits, 0.636 
and 0.646. We identify these limits with the random loose packing volume fraction <j> r \ p and the 
random close packing volume fraction <f> rcp of frictionless spheres, respectively. 
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INTRODUCTION 

Packing problems are among the most ancient. For in- 
stance, about 2200 years ago Archimedes faced the prob- 
lem of counting the number of grains of the beaches of its 
home town, Syracuse, and succeeded in demonstrating 
that this number is finite (Archimedes, The sand rock- 
ener). Another eminent Greek scientists, Apollonious of 
Perga, is also renowned for his works on disk and sphere 
packings. Modern works on disordered packing of grains 
dates back to the experiments conducted by Bernal[l| 
and Scott who prepared packings of monodisperse 
spheres using different protocols, and were able to mea- 
sure their volume fraction 0, defined as the fraction of 
the total volume occupied by the spheres. Their re- 
sults indicated the existence of an upper and of a lower 
bound for the volume fraction of stable disordered ar- 
rangements of spheres, named 'Random close packing', 
4> lcp , and 'Loose random packing', <p T \ p . Bernal estimated 
rcp = 0.63 ± 0.07 and <p lip = 0.60 ± 0.01, in the infinite 
system size limit. Regarding the existence of any first 
principle definition of these two bounds, Bernal[l[ specu- 
lated that 'The figure for the occupied volume of random 
close packing must be mathematically determinable, al- 
though so far we known undetermined'. Conversely, he 
questioned the existence of a first principle definition of 
the lower bound (f> T i p : 'The mathematical status of phys- 
ical random loose packing is not so evident'. 

Subsequent works clarified that (f> T \ p depends on the 
Coulomb friction coefficient. Onoda and Linigcr[3] op- 
eratively defined r i p as the smallest volume fraction at- 
tainable letting the particle sediment under gravity. This 
lower bound is obtained in the limit of zero sedimentation 
velocity, when particles fall in a very high viscous fluid. 
Onoda and Liniger[3j estimated r i p = 0.555 ± 0.005, 
while a more recent work[4] suggests r i p = 0.550±0.001. 

Since (f> r \ p depends on the friction coefficient, it is of 
interest to consider its value in the frictionless case. Pub- 
lished results seem to indicate that, in absence of friction, 
Consider, for instance, the 'jamming phase 
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diagram' introduced by Liu and Nagel[5j, illustrating the 
region of the temperature, volume fraction and stress 
space where jammed (mechanically stable) systems are 
found. Along the volume fraction axis, i.e. at zero tem- 
perature and zero applied stress, the transition between 
the unjammed and the jammed phase is marked to occur 
at a single value of the volume fraction, known as the J- 
point and later identified with the random close packing 
volume fraction (/> rcp . Since for <j) < rcp there are not 
jammed states, the diagram suggests that at zero fric- 
tion (/> r i p coincides with (/> rcp . The same conclusion could 
be drawn from the results of O'Hern and coworkers 16, 7], 
who numerically generated jammed packings of soft fric- 
tionless spheres using the conjugate-gradient protocol. 
In the hard sphere limit, they found jammed packings 
only at the volume fraction 4>* = 0.639 ± 0.01, identi- 
fied with (f> Icp . Being the only volume fraction at which 
jammed packings are found, it seems obvious to also 
identify <jf with r i p , concluding that at zero friction 

0rcp = 0rlp- 

In this manuscript we give evidence that, at zero fric- 
tion, </> r i p < rcp . An indication suggesting this possibil- 
ity comes from the comparison of numerical results found 
by different research groups, which have used slightly dif- 
ferent algorithms to prepare jammed packings of friction- 
less particles. For instance, using the conjugate gradient 
method [flQ, O'Hern et al. suggested 4> rcp ~ 0.639, while 
using a packing inflation algorithm Zhang and MakseQ 
obtained <f> lcp ~ 0.645 (with an error of the order of 
10~ 5 ). A very close value, 0.644, was also reported in 
previous works [sHU]]. These estimates are close, but not 
consistent within the reported errors. The discrepancies 
could be in principle attributed to one of the following 
causes: 

A: Finite-size effects. One or both estimates are 
wrong as affected by finite-size effects. 

B: Ordering. The upper bound for the volume fraction 
of grain packings is that of the FCC crystal, 4>fcg — 
0.74. One could therefore speculate that jammed 
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packings with volume fraction above 
contains some crystalline patches. 
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Disordered jammed packings occur 
in a whole volume fraction range, at least varying 
from 0.639 to 0.645. 

The possibility A) must be excluded, as finite-size are 
known to influence the jamming volume fraction leading 
to a smaller estimate of the critical packing fraction [ioj]. 
However, the smaller of the above estimates, </> rcp ~ 
0.639, has been obtained via a careful study of the in- 
finite system size limit 0, 0]. Here we introduce an al- 
gorithm able to generate jammed packings in a large 
volume fraction range (enclosing the range 0.639-0.644), 
and show that point B must also be excluded. In fact, 
altough the concept of random close packing of spheres 
is ill-defined [12[ due to the absence of precise definition 
of 'randomness', nevertheless our results indicate that it 
is possible to generate jammed disordered state with no 
cristalline patches up to the volume fraction <f) — 0.646, 
which is our estimation of the random close packing vol- 
ume fraction. Accordingly, we suggest that point C above 
is correct, i.e. that the jamming transition of frictionless 
spheres occurs along a whole volume fraction range. 

The manuscript is organized as follows. In Sec. „ we 
describe our numerical model and the protocol used to 
generated jammed packings of frictionless spheres. In 
Sec„ following Onoda and LinigerQ, we define as ran- 
dom loose volume fraction <p r \ p the volume fraction which 
is attained in the limit of infinitesimally slow energy min- 
imization, and determine its value for frictionless spheres 
in no gravity. We show that it is possible to generate 
un-jammed packings with a volume fraction which is 
higher than (j} T \ p in Sec. „ and describe their degree of 
order in Sec. „ This analysis clarifies that disordered un- 
jammed states cannot be obtained when the volume frac- 
tions overcomes a threshold, we identify with the random 
close packing volume fraction, (p lcp > 4> r \ p . Open ques- 
tions and future directions are presented in Sec. „ 




Generalized coordinate 



FIG. 1: Illustration of the possible states resulting from an en- 
ergy minimization protocol. These states can be jammed with 
a finite energy, jammed with zero energy (transition states), 
or unjammed with zero energy. 
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0. 



j/D) 2 if r < D, 
otherwise, 
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where r; specifies the position of particle i, rij = — rj\ 
is the distance between particles i and j. D, m and e are 
our units of length, mass and energy, respectively. The 
elastic E energy of this system is 



E({r,r}) 
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E 



v(nj), 
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where fi is the velocity of particle i, and N the number 
of particles. 

Jammed configurations correspond to energy minima 
of the system. Here we consider a configuration jammed 
when the mean energy per particle is higher than a small 
threshold, 10 _10 e. In principle, the definition of jamming 
we are considering allows for the presence of particles 
with no contacts (rattlers). However, in the infinite sys- 
tem size limit their concentration is known to vanish [7|. 



PROTOCOLS 

In this section, we introduce the numerical model, 
give the definition of zero pressure jammed packings 
of spheres, and describe the protocols used to generate 
jammed packings of spheres. 



Numerical model and zero pressure jammed packing 
of spheres 

We consider a system of monodisperse frictionless soft 
spheres of diameter D, and mass m, interacting with an 



Protocols 

Jammed packings of spheres correspond to minima of 
the energy of the system E = v ( r ij)- We sample 

these minima solving the equation of motion of the sys- 
tem in the presence of a viscous damping term, as if the 
particles were immersed in a viscous fluid, starting from 
a random (infinite temperature) configuration. Each par- 
ticle evolves according to: 

dE . 

™>n = -n ^i' ( 3 ) 

on 

where the parameter r\ plays the role of a viscosity. As 
time goes on, the total energy of the system decreases 
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because of the dissipative term, and the dynamics even- 
tually stops (the kinetic energy vanishes) . As illustrated 
in Fig. [TJ the resulting state will be either jammed with 
zero energy, or jammed with a finite energy. The figure 
also illustrates the existence of transition states, jammed 
states of zero energy. Energy minima may be found using 
other protocols, as for instance simulated annealing, in 
which the system is coupled to a heat bath whose tem- 
perature is decreased until it vanishes. Alternatively, one 
could consider the minimization of the energy of the sys- 
tem using as initial state an equilibrium configuration at 
a given value of the temperature: this is the protocol used 
to unveil the features of the energy landscape sampled by 
glass-forming liquids [2l|, [22 1. 

Qualitatively, the role of the parameter rj in the energy 
minimization procedure is easily understood. When 77 is 
high, the kinetic energy of the system is small, which im- 
plies that the system is less able to escape from the energy 
basins it visits. Accordingly, the larger r\ the smaller the 
region of the configurational space the system explores 
before getting trapped in a energy basin. In particular, 
in the limit r\ — > 00 the system jams in the first minimum 
it enters. Conversely, at a finite value of rj, the system 
may reach a minimum of the potential energy with a fi- 
nite value of the kinetic energy, and may therefore be 
able to escape from it overcoming the confining energy 
barriers. The average value of the energy of the minima 
reached with a given value of 77 is therefore expected to 
decrease as r\ decreases. 



RANDOM LOOSE PACKING 

Onoda and LinigerQ operatively defined r i p as the 
smallest volume fraction attained via sedimentation pro- 
tocols under gravity. Since sedimentation is an energy 
(gravitational energy) minimization protocol, this opera- 
tive definition can be extended to the absence of gravity, 
simply defining </> r i p as the smallest volume fraction of 
jammed states obtained via energy minimization proto- 
cols, the initial state being a random one. 

This jammed state of large volume fraction is obtained 
solving the equations of motion Eq.[3]in the limit 77 — > 00, 
in the infinite time limit. Operatively, this state can be 
also obtained via more tractable numerical procedures 
able to find energy minima, namely the steepest descend 
method, or the equivalent but computationally more per- 
formant conjugate gradient (CG) method. Here we use 
the Fletcher-Reeves CG al gori thm, as implemented in 
the GNU scientific library [15(. The CG is character- 
ized by two parameters. First, there is a tolerance small 
parameter, which is used to decide when the minimiza- 
tion along a given direction of the configurational space 
stops. Afterwards, the algorithm tries to minimize the 
energy moving the system along a conjugate direction 
of the phase space. In the minimization of the energy 
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FIG. 2: Probability Pjv of obtaining a jammed configuration 
minimizing the energy via the CG protocol, as a function of 
the volume fraction (j>, for different values of N. Each point 
is estimated performing 200 independent CG minimizations. 
For each value of N , the data are fitted by a scaled error 
function with inflection point (/>n and standard deviation ajv 
(plain lines). The inset shows the corresponding data collapse. 



E = v ( r ij) w e have not found a dependence of the 

results on this parameter, when this is small enough. 
Second, there is a threshold parameter Se: the algo- 
rithm stops when the energy difference between two suc- 
cessive iterations is smaller than this threshold. Here 
we set Se = 0, meaning that the algorithm stops when 
the energy variation in successive iterations is smaller 
than our numerical precision. We note that zero-energy 
jammed packings obtained via the CG method have been 
previously 0, 0| identified with </> rcp , while here we iden- 
tify them with </> r i p . We comment on this point later on. 

To determine r i p in the infinite system size limit we 
have performed a finite-size scaling, considering systems 
with a number of particles N varying from 170 to 4096. 
For each value of N, we considered different values of the 
volume fraction 0, performed 200 independent energy 
minimization CJ protocols, and measured the fraction 
P/v(0) of these minimization procedures yielding jammed 
configurations. To this end, we have considered a config- 
uration as jammed when its mean energy per particles is 
greater that 10~ 10 e, unjammed otherwise. At small vol- 
ume fraction, all of the minimization procedures result in 
un-jammed states, and Pn{<P) — 0, while conversely at 
high volume fraction all of them result in jammed con- 
figurations, and Pn(4>) — 1. The volume fraction depen- 
dence of Pn (4>) is well described by an error function 



Pn{4>) = \ 



1 + erf 



2a% 



(4) 



where an d crjv are estimated via a least square fit. 
Fig. [5] shows the raw data for different N, and their fits 
with Eq. [U 

As N grows, 4>n approaches an asymptotic value (j)^, 
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FIG. 3: Main panel: as the system size increases, the in- 
flex 4>n of the Pjv(0) data approaches the asymptotic value 
4>oo = 4>*lp as a power law (Eq. [5]| . r i p is well defined in the 
thermodynamic limit, as the width of the error function fits 
(jjv vanishes as N increases (inset). 



estimated [?] via a power law fit of the 
in Fig. [3] the power law fit 
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data. We show 



(5) 



where d = 3 is the dimensionality of the system, c^oo = 
0.636 ± 0.001 and v = 0.46 ± 0.06. The value of ^ is 
close to the one (0.639 ± 0.001) reported;/] by O'Hern 
et al., while the discrepancy in the estimation of v are 
more marked. Our estimation of </>oo is compatible with 
the smallest value of the volume fraction at which the 
pressure of monodisperse hard spheres has been recently 
found to di verge [19(. cf>oo is our best estimation of the 
random loose packing volume fraction in the infinite sys- 
tem size limit, <f> T \ v = <f>oo = 0.636 ±0.001. This value can 
be considered sharply defined, as the standard deviation 
fj 7v vanishes as a power law as N increases, as shown in 
Fig. GS (inset). 

We note that, within our numerical accuracy an expo- 
nential law describes the data dependence of 4*n on N 
equally well (in terms of the x _ square), an d provides a 
slightly different estimation of r i p , r i p ~ 0.6355. 



ABOVE <^i P 

We have determined the value of 4> v \ v finding the energy 
minima of the system via a procedure which is equivalent 
to the solution of the equation of motion (Eq. [3J , in the 
?/ — > oo limit, using as initial state a random one. Here 
we describe the results obtained when the minima of the 
system are obtained solving Eq. [3] in the presence of a 
finite value of the viscosity 77, until the dynamics halts. 
The initial state of this minimization procedure, however, 
is not a true random one. This is so because the correct 
simulation of the relaxation dynamics of random states, 



FIG. 4: The probability Pn of obtaining a jammed configu- 
ration minimizing the energy of the system via the numerical 
solution of Eq. as a function of the volume fraction (j>. The 
figure shows the results obtained for two different values of 
the viscosity 77, and for several values of the system size, as 
indicated. Each point is estimated performing 100 indepen- 
dent minimization procedures. Plain lines are fit to a scaled 
error function, Eq. |4] 



which may have very high elastic energy, requires the use 
of a very small numerical integration timestep, and is 
therefore too computationally expensive. We therefore 
follow Zhang and MakseQ, and use the following proto- 
col. We first prepare the system in a low volume fraction 
state with zero energy (no particle contacts). The size 
of the particles is then quickly inflated until the desired 
value of the volume fraction is reached (we increase <f> 
linearly in time). Afterwards, Eq.[3]is numerically solved 
until the dynamics halts. Accordingly, this procedure de- 
pends on two parameters, the rate T at which the size of 
the particles is varied, and the viscosity rj. We use the 
value of T considered by Zhang and MakseQ. 

For each value of the number of particles N, we have 
considered different values of the damping parameter 77, 
and different values of the volume fraction <fi. For each N, 
4> and 77 triple, we have performed 100 simulations of the 
relaxation process, and determined the fraction of these 
simulations which resulted to be jammed, -P/v(0, 77). We 
consider a system to be jammed when the elastic energy 
per particle is greater that 10~ 10 £. Raw data for two 
values of 77 are shown in Fig. 2J Fitting P/v(</>, ??) using 
an error function (Eq. [4j we have estimated the inflec- 
tion point 4>n(ii) an d the variance (7^(77). As N grows, 
4>n(i]) quickly approaches its asymptotic value, as shown 
in Fig. [5j We have therefore identified the N —> 00 limit 
i'ooi'n) with the value obtained with N = 8000 particles, 
the largest number of particles with have considered. 

The data of Fig. |4] clearly show that by changing the 
viscosity parameter 77 used in the minimization of the 
energy, different values of the volume fraction (j)^ are 
obtained. The dependence of </>oo(??) on 77 is shown in 
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FIG. 5: Crossover volume fraction 4>n as a function of the 
number of particles, for different values of the viscosity 77. 
From bottom to top, 77 = 20,10,8,5,2,1,0.1. (j>N quickly 
reaches an 77 dependent asymptotic value as N grows. 
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FIG. 6: The volume fraction (pea of the infinite system size 
limit as a function of the viscosity parameter r\ used in the en- 
ergy minimization procedure, in a log-linear (main panel) and 
in a linear scale (inset). At a given value of 77, if <j> < 4>oo(i]) the 
minimization procedure results in an unjammed state, while 
conversely a jammed state is obtained. Here, we approxi- 
mated with 0jv=8ooo- Errors are the standard deviations 
o"jv=8000 estimated via a least square fit of the jamming prob- 
ability Pn{4>,ti) with the error function of Eq.[4] 



Fig. [6l 0oo monotonously decreases as 77 increases, and 
the limit r) -> 00 yields (f>oo(T, 77 -> 00) = 0.634 ± 0.001, 
a value which is close to our estimation of the random 
loose packing volume fraction, </> r i p = 0.636 ±0.001. Rig- 
orously, one should expect </> r i p = </>oo(r — > 00,77 — > 00), 
as only in the r —> 00 the initial state of the minimiza- 
tion procedure is a random one. As 77 decreases, 0oo 
increases, and reaches values which are well above any 
past estimate of the random close packing volume frac- 
tion. The expected emergence of ordering in these high 
volume fraction states is described in the next section. 



RANDOM CLOSE PACKING 

The introduced energy minimization protocol allows 
the generation of unjammed packings with a volume frac- 
tion <j) which is well above any past estimation of the ran- 
dom close packing volume fraction, as shown in Fig. |6] 
Here we quantify the degree of order of these states con- 
sidering systems with N = 8000 particles, which Fig. [5] 
showed to be representative of the infinite system size 
limit. 

we quantify the de- 
whose defini- 



Following previous works 23 2£ 
gree of order focusing on the parameter Qi 
tion requires the introduction of bonds between particles. 
Here we consider two particles as bonded if they share a 
face of the Voronoi tessellation of the system. A bond 
between particles i and j defines a vector ry = (r, 9, 4>). 
Each bond can be therefore associated a whole set of 
spherical harmonics Yi m (6, (j)). By combining the values 
of the spherical harmonics associated to all bonds of a 
given particle, it is possible to associate to each particle 
a scalar parameter whose value de pen ds on the shape of 
its Voronoi cell. This parameter is|23j 



47T 

21 + 1 



1/2 



(6) 
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where the average (Yi m (9, (j))) is performed over all bonds 
of the particle. The sum over m assures the indepen- 
dence from the chose reference frame, i.e. makes Qi ro- 



tationally invariant. Previous works 23 25] have clarified 
that the most convenient value of I for the study of the 
emergence of crystallization in systems of hard spheres 
is I = 6, which is the lowest nonzero Qi in common with 
the icosahedral symmetry and cubic symmetry. Ordered 
structures are characterized by the frequent occurrence 
of typical local structures (i.e. particular shapes of the 
Voronoi cells), and are therefore characterized by a prob- 
ability distribution of P(Qe) with distinct peaks. For 
instance, all Voronoi cells of the simple cubic crystal are 
equal, and P{Qe) = 5(Q 6 - Qsc), with Q sc ^ 0.35. 
More complex crystals are characterized by the presence 
of Voronoi cells with few shapes, and lead to a proba- 
bility distribution P{Qe) with more peaks. For this rea- 
son investigating the evolution of P(Qe) is a convenient 
way to monitor the emergence of order into a system, as 
for instance the crystallization process of a supercooled 
liquid 25j. In particular, growing peaks at Qhcp — 0.48, 
Qbcc — 0.51, and Qfcc — 0.57 reveal the emergence 
of crystalline patches with the hexagonal close pack, the 
body centered cubic, and face centered cubic symmetry, 
respectively. It must be noted, however, that the as- 
sociation of a particular value of Qe to a precise local 
crystalline structure is not straightforward in the pres- 
ence of noise. For instance, the peak at Qq ~ 0.51 may 
also correspond to a FCC Voronoi cell with 13 faces 25 1. 
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FIG. 7: Probability distribution of the order parameter Qq 
of unjammed packings for <j> = 0.635, 0.646, 0.654, 0.670 and 
0.692. The curves are obtained at different values of rj, as 
indicated. As the volume fraction increases, peaks at values 
of Qe known to characterize different crystalline structures 
appear. 

We have investigated the probability P(Qe) of the un- 
jammed configurations generated minimizing the energy 
with the described protocol. The distribution depends 
both on the volume fraction <f>, and on the viscous param- 
eter rj. Representative distributions are shown in Fig. [7] 
At low volume fraction, P(Qq) has a smooth shape, while 
peaks are clearly present at high volume fraction. These 
peaks signal the occurrence of local arrangements typical 
of the HCP, of the BCC (or FCC with 13 faces) and of 
the FCC crystal. Other peaks correspond to the presence 
of many nearly-cristalline Vorono'i cells. 

From the probability distribution P(Qe) it is possible 
to extract a scalar order parameter which quantifies the 
degree of order by measuring the hight of the peaks with 
respect to the 'base' of the distribution. Here we suggest 
to first identify the n p values Q\ , k = 1 , . . . n p where 
the distribution P(Qq) has peaks, and then to define the 
order parameter as 

»p 

AP = ]T P(Ql) - - [P(Ql + SQ) + P(Ql - 6Q)} (7) 

where the value of SQ is irrelevant as long as this is larger 
than the width of the peaks, and smaller than the dis- 
tance between consecutive peaks. This order parameter 
compares the height of the peaks of the P{Q§) distribu- 
tion at values of Qq characterizing the ordered structures 
with estimates obtained via a local linear approaximation 
of the distribution. If the distribution has no peaks sig- 
naling the presence of ordered structures, then AP ~ 0, 
while AP > signal the presence of crystalline patches. 
We have computed AP fixing SQ = 0.1, and n p = 3, 
where Q\ = Qhcp, Q\ = Qbcc, Q3 = Qfcc- The inclu- 
siong of peaks observed at smaller Q§ (e.g Qe — 0.43 or 



FIG. 8: Order parameter AP of unjammed states as a func- 
tion of the volume fraction, for different values of the viscous 
parameter 77 = 0.01, 0.05, 0.6, 0.8, 1.0, 5.0 and 8.0 used in 
the minimization procedure. For <p < 0.646, there are values 
of rj leading to disordered unjammed states with AP ~ 0. 
Conversely, for (f> > 0.646 we find AP > for all values of 
Tj, suggesting that all unjammed states with (f> > 0.646 have 
crystalline patches. 

Q 6 ~ 0.38) is irrelevant, as these peaks correspond to 
slightly irregular crystalline Vorono'i cells, and are never 
observed alone. 

The dependence of AP on the volume fraction, for 
different values of 77, is shown Fig. and reveals the 
presence of a transition occurring at a volume fraction 
cj) ~ 0.646. Regardless of the value of 77, all unjammed 
states with <fi > 0.646 have some degree of order. Con- 
versely, when cj) < 0.646, high values of rj lead to disor- 
dered unjammed states. 

The jamming and the ordering properties of the sys- 
tem are summarizied in the <p-rj diagram shown in Fig. [9] 
This diagram clarifies that the jamming transition bew- 
teen disordered states occurs in a volume fraction range, 
with extrema r i p c_ 0.635 and </> rcp ~ 0.646. The jam- 
ming transition line also extends to values of cf> > rcp , 
as shown in Fig. [6l but at this high valued of the volume 
fraction ordering is found in the separated unjammed 
and/or jammed states. Fig. |9] clarifies that the tendency 
towards the formation of crystalline patches is higher 
the smaller the viscosity parameter. In this sense, the 
conjugate gradient protocol^, 0] is the lesser prone to- 
wards crystallization. The question [ilj] whereas there are 
jammed states below </> r i p , which may be obtained using 
different protocols, will be discussed in the next session. 

DISCUSSION 

We have introduced a protocol able to generate 
jammed zero pressure disordered packings of frictionless 
spheres with a volume fraction varying in a whole range. 
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FIG. 9: Jamming and ordering properties as a function of 
the volume fraction and of the viscosity parameter. The 
shaded area covers the region of the (p-r/ plane where crys- 
talline patches are found in the generated packings. The large 
full circles identify the jamming transition line between fluid 
like disordered states and disordered solid states. The ex- 
trema of this line are the random loose (77 — > 00) and the 
random close packing volume fraction, respectively. 



The obtained volume fraction depends on the value of a 
viscosity parameter 77. The lower extremum of this vol- 
ume fraction range is obtained solving the equations of 
motion of the system in the quasistatic limit 77 — > 00, 
starting from a random initial condition. Following On- 
oda and Linigerjsj, the quasistatic minimization of the 
energy can be considered as the operative definition of 
the random loose volume fraction. Our results suggest 
(f) r i p ~ 0.636. 

As the viscosity parameter 77 decreases, the volume 
fraction of the jammed configurations increases. Intro- 
ducing an order parameter based on the probability dis- 
tribution of finding Vorono'i cells with peculiar shapes, we 
show that our numerical protocols generate un-jammed 
disordered states up to a volume fraction rcp ~ 0.645, 
we identified with the random close packing volume frac- 
tion. We cannot exclude that its value may change if one 
considers a different order parameter. For </> > lcp all 
unjammed states appear to contain crystalline patches. 
Zero pressure disordered jammed frictionless spheres can 
therefore be found in a whole volume fraction range, at 
least extending from </> r i p to </> rcp . This scenario is consis- 
tent with recent numerical results investigating the jam- 
ming transition of thermal systems 
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An insight on the origin of these limiting values of the 
volume fraction is obtained considering how many differ- 
ent jammed zero pressure states exists at any given value 
of the volume fraction 4>. This leads to the introduc- 
tion of the density of disordered jammed zero-pressure 
states 13, [l4 1, fi(</>), or of the entropy S{<p) = logf2(0). 
We have recently investigated this quantity in a two di- 
mensional model(l3l|. finding an non- monotonous entropy 



FIG. 10: Schematic [13|] dependence of the entropy S(4>) = 
log £l(4>) on the volume fraction, where Q(<j>) is the density of 
disordered jammed states. 



S((j}) which is zero at small 0, then increases with 4> up to 
its maximum value, and finally decreases with <p until it 
vanishes at high <fi, as scheamtically illustrated in Fig.HUl 

When minimizing the energy in the infinite vicosity 
limit starting from random initial configurations, it is 
reasonable to assume as a first approximation that one 
finds all jammed states with the same probability (more 
precisely, one should consider that the probability of find- 
ing a minima is proportional to the width of its energy 
basin). If this is so, then the volume fraction obtained in 
the 77 — > 00 limit is the one where the entropy S((/)) has 
a maximum, which should be therefore identified with 
the random loose packing volume fraction. Likewise, 
since our results suggest that there are not disordered 
unjammed packing with (j) > </> r cp, then </> rcp could be 
associated with the value of the volume fraction where 
the entropy vanishes, at high volume fraction. Accord- 
ingly, the non-monotonous variation of the entropy with 
the volume fraction 13[ indicates that there exist dis- 
ordered zero pressure jammed states with volume frac- 
tion <fi < r i p . The smallest volume fraction of these 
states should correspond to the random very loose vol- 
ume fraction [l^l, rv i P - 

Open questions ahead include the developing of proto- 
cols to generate disordered jammed states of with volume 
fraction (f> < </> r i p , as well as the clarification of the vol- 
ume fraction dependence of mechanical and geometrical 
properties of jammed packings. 
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